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Abstract 

The asymptotic state discrimination problem with simple hypotheses is considered 
for a cubic lattice of bosons. A complete solution is provided for the problems of 
the Chernoff and the Hoeffding bounds and Stein's lemma in the case when both 
hypotheses are gauge-invariant Gaussian states with translation-invariant quasi-free 
parts. 

1 Introduction 

Assume that we know a priori that the state of an infinite lattice system is either pi [null 
hypothesis H ) or p 2 (alternative hypothesis Hi), and we want to decide between these two 
options, based on the outcome of a binary measurement on a finite part of the system. 
Obviously, there are two ways to make an erroneous decision: to accept Hq when it is false 
(error of the first kind) and to reject it when it is true (error of the second kind). In general, 
one cannot make the corresponding error probabilities to vanish, but they are expected to 
vanish in the limit as we increase the size of the local system on which the measurement is 
made. 

Hypothesis testing results [H [21 El HH C21 E31 [IH HSl HSl UHl EI] show that, in various 
settings, the optimal error probabilities actually decay exponentially, and the exponent of 
the optimal decay rate can be expressed as a certain generalized distance of the states p\ 
and p2, depending on the concrete setting of the problem. The most studied cases are the 
problems of the Chernoff and the Hoeffding bounds and that of Stein's lemma, and the 
corresponding generalized distances are the Chernoff and the Hoeffding distances and the 
relative entropy, respectively. Apart from giving computable closed expressions for the error 
exponents, the importance of these results lies in providing an operational interpretation 
for the corresponding generalized distances, which in turn yield alternative and heuristically 
very transparent proofs for their monotonicity under stochastic operations [U [18] . 



1 Electronic mail: milan.mosonyi@gmail.com 



1 



The first such result in the quantum setting was obtained by Hiai and Petz [13] and 
completed later by Ogawa and Nagaoka [21], solving the problem of Stein's lemma in an 
i.i.d. setting (i.e., when p\ = and P2 = (p^)^ 00 are translation-invariant product 

states) on a one-dimensional spin lattice. This result was later extended to certain correlated 
situations as well as to higher dimensional lattices [3l[T5]. The recent findings of Nussbaum 
and Szkola [19] and Audenaert et al. pQ created renewed interest in hypothesis testing 
problems, and their methods were successfully applied to solve such problems in various 

settings [DElEldllinillSlinillBlIIH]. 

The study of the Chernoff bound for identical copies of one-mode Gaussian states in 
bosonic systems was initiated in [6] , where an explicit formula for the Renyi relative entropies 
was provided, which was later generalized for n-mode states in [23]. Stein's lemma for 
identical copies of one-mode gauge-invariant Gaussian states was treated in [10]. Here we will 
study the hypothesis testing problem for gauge-invariant Gaussian states with translation- 
invariant quasi-free parts on an infinite bosonic lattice, and we give a complete solution for 
the problems of the Chernoff and the Hoeffding bounds and Stein's lemma in this setting. 

The structure of the paper is as follows. In Section [27T1 we give a more technical introduc- 
tion into hypothesis testing and in Section 12.21 we overview the basic facts about Gaussian 
states that we will use in the rest of the paper. In Section 13.11 we prove the existence of 
various asymptotic quantites, including the mean Chernoff and Hoeffding distances and the 
mean relative entropy, and in Section 13.21 we show that these quantities give the optimal 
decay rate of the error probabilities in the corresponding settings. In computations with 
Gaussian states, we will need some basic facts about Fock operators that we collect in a 
separate Appendix. 

2 Preliminaries 

2.1 Hypothesis testing on infinite lattice systems 

Consider a lattice system on a //-dimensional cubic lattice Z u . We assume that the observ- 
ables of the system span a C*-algebra A and the shift operations on the physical lattice lift 
to an automorphism group 7k, k 6 Z" on A, such that the observable algebra A\ corre- 
sponding to a finite A C Z" is generated by {7k (*A{o}) , k G A}. Typical examples are spin 
lattices, when A = £g>kez"£> {^ji an< ^ fermionic/bosonic lattices, when A is the CAR/CCR 
algebra on the anti-symmetric/symmetric Fock space on I 2 {~L V ). We also assume that for all 
finite Ac Z" there exists a Hilbert space TC\ such that TC\ = ®keA^o an d -4a C B(H\), 
which is satisfied by all the examples mentioned above. States of the infinite lattice system 
are described by positive linear functionals on A that take the value 1 on the unit of A. 

Assume that we know a priori that the state of the infinite lattice is either p\ (null 
hypothesis Hq) or P2 (alternative hypothesis Hi). To decide between these two hypotheses, 
we can make binary measurements on finite parts of the system, and for simplicity we assume 
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these finite parts to be ^-dimensional cubes C n := {k £ 1_ v : k±, . . . , k v — 0, . . . , n — 1}. A 
test on C n is an operator T £ ^4c„) < T < 7 n , that determines the binary measurement 
with measurement operators T and I — T. If the outcome corresponding to T occurs then 
H is accepted, otherwise it is rejected. The error probabilities of the first and the second 
kind are then expressed as 



a 



l (T):=p 1 (I-T) and (3 n (T) := p 2 (T). 



As noted in the Introduction, the error probabilities are expected to decay exponentially 
if we let n go to infinity and choose the measurements in an optimal way. Here we will be 
interested in the exponents 

dpi II P2) ■= sup 1 liminf — - log (a n {T n ) + (3 n (T n )) 1 , 

{Tn} [ rwoo n v J 

c(pi||p 2 ) := supilimsup — -log (a n (T n ) + (3 n {T n )) \ , 

{T n } I n^oo U ) 

c(pi 1 1 p 2 ) := sup { lim — - log (a n {T n ) + (3 n (T n )) 1 , 

{Tn} [n^co n v J 

corresponding to the problem of the Chernoff bound, 

M r IPillP2) := supiliminf — \\og(3 n {T n ) limsup — loga n (T n ) < -r 1 , r > 0, 



h(r\ Pi 1 1 P2) := sup \ lim sup log /3 n (T n 

{Tn} I n ^ n v 



limsup — loga n (T„) < -r > , r > 0, 



M r IPiHP2) : = sup { lim — ^-log/5 n (T„ 

{T„} n," 



limsup — loga n (T n ) < — r > , r > 



n— >oo 



corresponding to the problem of the Hoeffding bound, and 

s(pi||p 2 ) := supiliminf — -log/3 re (T n ) lima n (T„) = ol 

{T n } I n ^°° ^ J 



{^n} 



s(Pi 1 1 P2) := sup <^ lim sup log /3„(T, 



s(Pi||P2) := sup I lim — \ \og(3 n {T n 

{Tn} [n^oo n v 



lim a n (T n ) = 
lim a n (T n ) = \ , 



n^oo 



corresponding to Stein's lemma, respectively. The suprema are taken with respect to se- 
quences of tests, with T n £ Ac„, n £ N. Obviously, c(pi || p 2 ) < c(Pi \ \ Pi) < c(pi || p 2 ), 
fr(r| pi || p 2 ) < h(r\ pi || p 2 ) < h(r\ p 1 \ \ p 2 ), r > 0, and s(p 1 || p 2 ) < s(pi || p 2 ) < s(pi 1 1 p 2 ). 
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Assume for the rest that the restrictions p^ onto Ac n are given by density operators p^ 
on Hc n , i-e., there exist trace-class operators p^ 1 on 7ic n such that p% (a) = Trpj^a, a G 
Ac n - In this case, 

a n (T) = TrpS n) (/ n - T) and (3 n (T) = Trp^T. 

If the supports of p^ and p^ are orthogonal to each other then pf^ and p^ can be 
distinguished perfectly in the sense that there exists a test T for which a n (T) = (3 n (T) = 0, 
and hence the hypothesis testing problem becomes trivial. To exclude this case, we will 
assume that the supports of p^ and p^ are not orthogonal to each other for any ?i6N. 

Due to Holder's inequality, the operators (p^)* (p^) 1 * are trace-class for each t G [0, 1], 
and hence the functions 

^W^logTr^)*^) 1 "', tG[0,l] (1) 

are well-defined for each n G N, such that — oo < ip n {t) < 0, t G [0, 1], n G N. It is easy 
to see that ip n is convex on [0, 1] for all n G N. All along the paper we use the convention 
0* := 0, t G R, i.e., we take powers only on the support of p£ ' . In particular, (pj^)° denotes 
the support projection of p£ . For each n G N, the Chernoff disance of and p^ is 
defined by 

cip^Wpt^^-^Mt). 

For any r > 0, the Hoeffding distance of p^ and p^ with parameter r is 

U/ | (n) I, (nK ~*T - ^ n (t) 

H(r\p\ '\\p\') := sup . 

o<t<i i — t 

For the rest we assume that the additional condition 

(pi )*' ^ = 1,2, are trace-class for all t G (0, 1] and n G N (2) 

holds. Then, it is not too difficult to see (by using the eigen-decompositions and Lebesgue's 
dominated convergence theorem) that ip n is continuous on [0, 1] and different iable in (0, 1) 
for each n G N. Moreover, if suppp^ < suppp^ then ip n (l) = 0, and 

H(0\ p^llP?) = 5-^(1) = lim<(f) = TrpS n) (logpf - ]og$>) =: S(p^ ||p^), 



t/i 

where d~ip n (l) is the left derivative of ip n at 1, and 5'(p^ n ' ) 1 1 p^) is the relative entropy of p^ 
and p^. Though assumption (T2J) is quite restrictive in general, it is automatically satisfied 
when the local Hilbert spaces are finite-dimensional (which is the case for spin lattices and 
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fermionic lattices) and also when pi and p 2 are Gaussian states of bosonic lattices, as we 
will see later. 

One can easily see that if the sequence of functions -^pipn converges uniformly to some 
function if; on [0, 1] then the mean Chernoff distance and the mean Hoeffding distances exist, 
and 

C M (pi\\ P 2) := ]im^C(p { r ) \\p^) = -minm, (3) 

n— >oo n 0<t<l 

H m {t\ pi II pa) := hm —H{n»r\ pj n) \\ p { 2 n) ) = sup ~ tr ~ ^ r > . (4) 

n->oo n v 0<t<l 1 — t 

If, moreover, suppp-j 71 < suppp^ , n G N, and lim n -^d~il) n (l) = d~ip(l) then the mean 

relative entropy Sm(pi\ \ P2) '■= lim„^oo ■^S{p ( f'^ \ \ p 2 ) and the mean Hoeffding distance 
with parameter exist, and 

S M {px\\p%) =#M(0|pi||p 2 ) = 0-^(1). (5) 

A complete solution to the problems of the Chernoff bound, the Hoeffding bound(s) and 
to Stein's lemma is obtained if one can show that the relations ([3]), (jlj) and ([5]) hold, and 

C(pi || p 2 ) = C(pi || p 2 ) = C~(pi I) p 2 ) = C M (pl\\p2), 

h( r \pi\\p2) = h(r\p 1 \\p 2 ) = h(r\p 1 \\p 2 ) = H M (r\ p 1 \ \ p 2 ) , r>0 
s(pi || pi) = s(pi (I p 2 ) = s(pi (I p 2 ) = S M (Pl\\P2)- 

This was done, for instance, for i.i.d. states on a spin chain [TJ [HI [TU HE! EHJ El], for quasi- 
free states on a fermionic lattice jTTj and for Gibbs states of translation- invariant finite-range 
interactions on a spin chain [121 apart from the identity Sm(pi \ \ Pi) = d~~ip{l) (which, 
however, seems to follow from the results of [20]). Partial results were also obtained for 
finitely correlated states on spin chains in [T2l [T5] . Stein's lemma was also proven for the 
case when pi is an ergodic state and p 2 is a translation-invariant product state on a spin 
lattice [3]. 

2.2 Gaussian states on the CCR algebra 

Let (H,a) be a symplectic space, i.e., if is a real vector space and a is a non-degenerate 
antisymmetric bilinear form (a symplectic form) on H, and let k be a positive real number. 
We say that a map W : H — > A to a C*-algebra .4 is a realization of the (k, a)-canonical 
commutation relations (CCRs) if A is generated by (H^(x) : x G -ff}, and 

W(x)* = W(-x), W(x)W(y) = e- iKa{x ' y) W(x + y), x,y G H. 

Obviously, kg is again a symplectic form, and hence the introduction of k may seem super- 
fluous in the definition. However, we follow this terminology in order to be as compatible 
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as possible with the various conventions appearing in the literature. For a more detailed 
treatment of the following, we refer the reader to (HI E2] • 

By Slawny's theorem, any two realizations W\ : H — > A\ and W2 '■ H — > A2 of the 
(k, o")-CCRs are isomorphic to each other in the sense that there exists a C*-algebra iso- 
morphism a : A\ — > A2 such that a o W\ = W%. Moreover, if the symplectic spaces 
(H\, «i(Ti) and (H2, ^2^2) are isomorphic to each other then any two representations of the 
corresponding CCRs are also isomorphic to each other. As a consequence, realizations of 
the (k, ct)-CCRs with the same a and different k's are all isomorphic to each other, and we 
denote by CCR(H, a) the C*-algebra generated by any such realization. Also, since any two 
finite-dimensional symplectic spaces of the same dimension are isomorphic to each other, 
so are the realizations of the corresponding CCRs. Hence, if H is finite dimensional with 
dimiJ = 2c? then one can assume without loss of generality that H = R 2d and a is its 
standard symplectic form 

d 

cr ((xi, . . . , x 2d ), (yi, • • • , V2d)) ■= ^2 ( x Mk+d - x k+d y k ) , 

k=l 

which is the usual choice in physical applications, and the parameter k is usually taken to 
be 1/2 or l/(2h). 

Here we will consider the situation when H = Ti, for some complex Hilbert space 
(considered with its real vector space structure) and a is its standard symplectic form 
^n( x iy) '■= Im(x, y) , x,y G H. Let V m H denote the mth antisymmetric tensor power 
of 7i, with V°7i := C, and let T{H) := ®„ =0 V m TC be the symmetric Fock space. For each 
x G 7i let xp := Ylm=o ~^ x ^ m e ^CH) denote the corresponding Fock vector (also called 
coherent vector or exponential vector). The Fock vectors are linearly independent and their 
linear span is dense in T(1~C). The Weyl unitaries W K (x), x G Ti. on Ti^K) are defined by 

W K (x)y F := e-HNI 2 -^,^ + ^ x)p ^ y e ^ 

and they are easily seen to give a realization of the (k, cr^)-CCRs. We denote the generated 
C*-algebra by CCR(ft). 

A state p of CCR(7Y) is a positive linear functional p : CCR(7^) — > C, that takes 
the value 1 on the unit of CCR(7i). The characteristic function of a state p is W R [p] : 
H — > C, H^fp^x) := p(W K (x)) , x G H. For any real inner product (i.e., positive definite 
symmetric real bilinear form) a satisfying 

a(x, y) 2 < a(x, x)a(y, y), x,y E H, (6) 

there exists a unique state p a on CCR(7i) with characteristic function 

W K [p a ](x) = p a (W K (x)) = e~^ x \ xeH. 
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Such states are called quasi-free. Note that is equivalent to the kernel (x, y) i— > ?/) + 
io~(x, y) being positive semidefinite. Obviously, 



PaM := p a (W K (y)*aW K (y)) , a G CCR(H) 

is again a state, with characteristic function 

W K [p a>y ](x) = e M^(».*)-f ? x ETi. 

States of this form are called Gaussian, and we will refer to 2/ as the displacement vector. 

The gauge group of CCR(7i) is the group of quasi-free automorphisms 7a, A G T := 
{z G C : |z| = 1}, defined by 7a (W k (x)) := W K (Xx), x G TL. A state p is gauge-invariant 
if p o 7 A = p, A G T. A Gaussian state p a)3/ is gauge-invariant if and only if a is gauge- 
invariant, i.e., a(Xx, Xy) = a(x,y), x, y G TL, X G T. Moroever, if TC is finite-dimensional 
then a Gaussian state p atV is gauge-invariant if and only if there exists a complex linear 
operator A > I such that 

a(x, y) = Re (Ax, y) , x, y G TL. (7) 

The operator A is called the symbol of p aiV . In the general case we say that p a has a 
symbol if there exists a complex linear operator A > / such that (JTj) holds. Note that 
if 7i is infinite-dimensional then having a symbol is a possibly stronger assumption than 
gauge-invariance . 

If p a>y is gauge-invariant and TC is finite dimensional then p a> y has a density operator, i.e., 
there exists a trace-class operator p a>y on JF(H) such that p a ,y( a ) — Tr p a>3/ a, a G CCR(7i). 
Moreover, the density operator can be expressed in terms of the symbol and the displacement 
vector in the form 

2dimW / j\_ j\ 

p a , y = W K (y)p a W K (yy, P*= det{l + A) {aTi) f > 



where Xp denotes the Fock operator corresponding to an operator X (see Appendix A ). For 
a proof, see e.g. [EJ Corollary 3.2]. Note that the eigenvalues of A coincide with the sym- 
plectic eigenvalues of a, and formula (1321) gives essentially the normal mode decomposition 
of the state in the above formalism. 

Consider now a bosonic lattice system on the //-dimensional cubic lattice Z u , such that 
to each physical site there corresponds one mode of the system. That is, the one-particle 
Hilbert space of the system is H := I 2 (Z u ) and its observable algebra is CCR(7i), the 
C*-algebra generated by the Weyl unitaries on J-(Tl). Let {l{k} : k G Z"} denote the 
standard basis of I 2 (Z y ), let TL n := span{l{ k } : k±, . . . , k v = 0, . . . , n — 1} and let P n be 
the projection onto H, n . For a bounded operator A on 7i, let A^ denote P n AP n , when 
considered as an operator on 7i n . The Hilbert space of the subsystem corresponding to a 
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cube C n := {k : k±, . . . , k v — 0, . . . , n — 1} is JF(7i n ) and its observable algebra Ac„ is the 
C*-algebra generated by {W K {x) : x G H n }, that is, CCR(7i n ). 

The shift operators are given by : l{j} i— > l{j+k}? j, k G Z u , and they induce the 
translation automorphisms of CCR(7Y), given by 7k : W K (x) i— >• W K (Skx) , i 6 H, k e Z". 
A state p on CCR(7i) is translation-invariant if p o 7 k = p, k G Z". Assume that p a is a 
quasi-free state given by a symbol A. Then, p a is translation-invariant if and only if A is 
translation- invariant, i.e., S^AS^ 1 = A, k G Z". 

Translation-invariant operators are diagonalized by the Fourier transformation 

T: l 2 (T) ^L 2 ([0,2tt)"), 0Flj)(x) := e i ^=i^ m> x G [0, 2tt)", j G Z 1 '. 

That is, A is translation-invariant if and only if there exists a bounded measurable function 
a : [0, 2n) u — > C such that A = J r ~ l M a J r , where M a denotes the multiplication operator by 
a. Let denote the convex hull of the spectra of A. We will make use of the following 

multivariate extension of Szego's theorem [9], that was proven in [IT] : 

2.1 Lemma. Let ai, . . . , a r be bounded measurable functions on [0, 2ti) u with corresponding 
shift-invariant operators Ai, . . . , A r . Then, 

hm i-TrA^) • . . . • / r (4»)) = / /i (a x (x)) •...•/, (a r (x)) rfx (8) 

for any choice of polynomials j\, . . . , f r . If all a k are real- valued then ([8]) holds when is 
a continuous function on £(^4^) for all 1 < k < r. In this case, the convergence is uniform 
on norm-bounded subsets of Ylk=i ^ (^(■^■fc))> where C (E(A fc )) denotes the vector space of 
continuous functions on E(Afc), equipped with the supremum norm. 



3 Hypothesis testing for Gaussian states 

Consider now the hypothesis testing problem described in Section 12.11 We will assume 
that pi = p ai ,yi and p 2 = p a2)V2 are both gauge-invariant Gaussian states, and, moreover, 
that their quasi-free parts p ak are translation-invariant and are given by the symbols A k = 
J r ~ 1 M ak J-', k = 1,2. Here, a k : [0,271")" — > [l,+oo) are bounded measurable functions. 

The restrictions p^,y k °f Pa k ,y k onto CCR(7i n ) are again Gaussian states, with symbols 
A k n ^ = P n A k P n and displacements y k n ^ = P n Vk- Moreover, phl, yh are given by the density 
operators pSl,y k = W K (y k n ^) p£} W K (y^)*, where p^ are the densities of the quasi- free parts, 
and 



2 " / A k n - 1 \ 1 / Qk ) 

- det (/ + 4"») U W + 'I f det (/ + Q<">) \qP + I ) ( 



8 



where 

Q k :=(A k -l)/2, R Kn -.= ^ an d N k>n := 1/ det (/ + Q k n) ) . (9) 

Qk + 1 

Note that Q k = J r ~ 1 M qk J-', where q k = {a k — l)/2 are non-negative bounded measurable 
functions on [0, 2ir) u , and p ak are uniquely determined by either of the functions a k , q k and 
r k '■= Qk/0- + Qk)- F° r later use, we define 

W n , t := R'Hr^R'H and w t (x) := r 1 (x)V 2 (x) 1 ^, t 6 R. 

Note that i4, n ^ R^ := (j^M^^ and W n , t ^ W t (n) := (J 7 ' 1 M Wt T) {n) in general. 
An easy computation yields that for < t < 1, 

[&J = M ^tP ft{a i^ y , M k , n , t : = 2*/ det [M4°)L 

where 

V( S ):=( S + 1)*±( S -1)*, / t ( s ) := ^±^l> *e(0,l), *>1, (10) 

and 

(z,y) := Re (f t (A^)x,y), x,yeH. 

That is, the powers {pr^ k yk \f are again the densities of Gaussian states up to a normalization 
constant, and hence condition (j2j) is satisfied. 

For some of the statements we will also have to assume that qi and g 2 are strictly 
positive in the sense that there exists some rj > such that ^(x) > r\ for almost every x, 
or equivalently, that Q k > This assumption ensures that the local restrictions p^, Vk are 
faithful for each n G N, or, in more physical terms, that the vacuum state does not appear in 
the normal mode decomposition of p& k for any n. Note that this notion of strict positivity 
is stronger then requiring q k (x) > for almost every x. 

3.1 Asymptotic distances 

Let A k = Jk,nj\ e k,nj) ( e k,nj\ be eigen-decompositions of the symbols, and assume that 
the eigenvalues jk,n,j are ordered so that Jk,n,j > 1 f° r 3 < r(k,n) and jk,n,j — 1 f° r 
j > r(k,n). By ©, this gives the eigen-decompositions R k . n = YJjtf* ^k,n,j\e k , n ,j) (ek,n,j\, 
where Xk,nj — 7 |?"' J +i > 0, j = 1, . . . , r(k, n). With the notations of (}33l) . we get the 
eigen-decompositions of the densities as 

oo 

Pctk,y = Xk,n,mjW K (y)e ktntm ){W K (y)e k:nt!R \, X k ,n,m '■= ^fe,n^fe,n,m- 

m=0 mi+...+m r =m 
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Following [15) , we define 

Pl,n {Wl,Wl) '■= Al, n ,m| (ei,n,m, e 2 ,n,m') | 2 , £>2,n (ZZl, ZZl') : = A 2 ,n,m' | (ei,n,m, e 2 ,n,m') | 2 , (lR,mf) G J„ 

where J n := {(m,m') G N r(1 ' n ) x N r ( 2,n ) : (ei jn>m , e 2j „ )ffi ) 7^ 0}. Then, p fcjn are positive 
measures on J n with Pk, n (Jn) < 1, and 

^ (Pai'.Vl) (Pwzm) = ^2 ^ln,m^2,n*m'l ( e l,ri,m, e 2 ,n,m') | 2 

(m,m')6Jn 

= ^2 Pi^im,™) 1 P2,n(ni,Ek) 1 ~ t , te[0,l]. 

(jn,m')eJ n 

Define 

ip n (t) := log ^ p 1>n (m,m') i p 2jn (m,^') 1- * , i 6 R, 

(m,m')eJ n 

using the convention log +00 := +00. By Holder's inequality, ip n is convex on IR for each 

Our first goal is to express ip n in terms of the symbols and the displacement vectors. To 
this end, let 

c n>t := exp ( - 2k{ [f t (A^) + f x _ t (A^)} ~ \ y» , t G (0, 1) (11) 

with y := y 2 — y\ and /j given in ffTUj) . Note that f t (s) > 1 for all t G (0, 1), s > 1 and 
lim^o /t(s) = +00 if s > 1. If A x 7^ / then 7^ J n , n G N, and we have lim^ [f t (A.[ n ') + 
/i_t (A 2 n ^)] = 0. If A\ = I (which is the case if and only if p ai is the vacuum state) then 
A^ = I n , n G N, and lim 4 ^ f t {A^} + fi- t (A^) — I n + A^\ By similar considerations 
with v4 2 and lim^i, we get 



t\o 1 exp ( - 2«(j/, + Ij y» , if Ai = J, 



if A 2 ^ /, 

vr n,£ \exp (2^<y, [Af ] + if A 2 = /, 

for all n G N. 

3.1 Lemma. For each t G [0, 1], 

(t) = log c M - 1 Tr log (Qi n) + I n ) - ( 1 - 1) Tr log (Q< n) + /„) - Tr log (l n -W n , t ). (14) 



10 



Proof. The Parseval formula for the Weyl transform [TB] tells that for any two Hilbert- 
Schmidt operators T 1; T 2 on FiTin), 



Tr T*T 2 =(-)*[ Tr(iy K (x)T!) Tr(W K (x)T 2 ) d\(x) 



'T~tn 

where A is the Haar-measure on 7i n , normalized so that cubes spanned by symplectic bases 
have measure 1. The choice 7\ := (p|™| yi \Y, T 2 := j/ 2 )) ' with t G (0, 1) yields, after 
some computation, 

2 d c n M l „ t M 2 „ t 



det [/ t (4 n) )+/i-t(4 B) )] 



(15) 



2 d c 



n.t 



det [(4^ + /n)'(4 W) + In) 1 '' - (4^ " 4)*(4 n) " In) 1 ^] 

C n ,t 



det [(g^ + iJ'^ + Jn) 1 -*- (qPY(qP) H 

From the last formula the assertion follows for t G (0, 1), and the cases t = and t = 1 can 
be verified by a direct calculation. □ 

3.2 Remark. Note that the above result gives that for all t G [0, 1], 

Mt) = logTr (p^J (pgj 1 -* = lo gCnii + logTr (/^J) 1 "* , 

i.e., the effect of the displacements only appears in the term logc nj f. 

The core of the above Lemma, formula (fT5|) was derived in [23] in the general (not 
necessarily gauge- invariant) case. If the two states have the same displacement then the 
above result might be strengthened and the proof reduces to a straightforward computation, 
as shown below: 

3.3 Lemma. Assume that y\ = y 2 =: y. For all t G IR such that W n> t < I n , 

Mt) = -tTrlog (Q^ + I n ) -(1-t) Trlog (Q< n) + /„) - Tr log (/ n - W n , t ). 

Proof. With the convention 0* := 0, t G R, the powers (paj!,yY are well-defined positive (not 
necessarily bounded) operators for any t G R, and 



By (1341) . the above operator is bounded and trace-class if and only if W n ± < I. 
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Let us define the trace of a positive operator to be +00 whenever it is not trace-class. 
With this convention, 



(Pi,n )\p 2 , n (m,rnf) )^ = Tr (fe] y ) t/2 (fetf-* (feW* = K> N £ Tr ( W ^)f , 

(m,m')eJ n 

from which the assertion follows. □ 

3.4 Remark. One can easily see that if q\ and q 2 are strictly positive then there exists some 
5 > such that W n ,t < I n for all £ G (—5, 1 + 5). In this sense, Lemma 13731 is an extension 
of Lemma [3.11 in the case when y\ = y 2 . 

Our next goal is to prove that the limit 

^(t):=lim^Vn(*) (16) 

exists for all t G [0, 1]. For this, the following simple Lemma will be useful: 

3.5 Lemma. 

M(q 1 ,q 2 ) := sup{||W n , t || , ||w t || : t G [0,1], n G N} < 1. 
Proof. Note that ||Qi n) || = \\P n QkPn\\ < \\Qk\\ = Ik/fclL, and hence, ||-R fc ,„|| < H^ILAl + 

IkfclL) = IkfclL- Thus for each t e [0) !]) 

\\W n , t \\ < \\Ri, n \f ||^2,n|| 1_ * < Ikill^ ||r 2 f < maxjllriH^ , IHI^} < 1, 
and the same bound holds for HtftH^, from which the assertion follows. □ 

3.6 Lemma. The sequence -^pipn converges uniformly on [0, 1] to 

m = -7A17 / log [(1 + fc(x))* (1 + ^(x)) 1 "' - ( gi (x)) 4 (^(x)) 1 "*] Ac. (17) 

If qi and q 2 are strictly positive then tp is twice different iable in (0, 1), and ip"(t) > 0, t G 
(0, 1) unless 51 (x) = q 2 (x.) for almost every x. 

Proof. By lemma 13. 1[ 

Mt) = log Cn.t - t Tr log (4 + gS n) ) - (1 - 1) Tr log (/„ + Q ( 2 n) ) - Tr log (J n - W n>t ) . 
Formulas ([II]), (H2D and ([13]) show that e"^ 2 < c n , t < 1, < t G 1, and hence, 

lim^-logc n . t = 0, te[0,l]. 

n n u 
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By Szego's theorem, 

1 1 f 

lim — Tr log (j n + Q^ ] ) = ——^ / log (1 + g fc (x)) dx. 



'[0,271-)" 

By lemma < W n>t < M(q 1 , q 2 )I n and < w t (x) < M(q u q 2 ) for all t G [0, 1], n G N, 
and almost every x. Consider the power series expansion log(l — x) = — Xlm=o which 1S 
absolutely and uniformly convergent on [0, M(q±, q 2 )}, and define Pn(x) : = — J2m=o ~i N G 
N. Then, 

< \\l0g (I n -W n , t )-p N (W n , t )\\ 



— v Tr log (/ n - W njt ) - — Tr p N {W n , t ) 



< max {|log(l — x) — Pn(x)\} ► 0. 

~ a-6[0,Ar(gi,ga)] W^oo 



Hence, it is enough to show that 

1 1 /" 

lim — Tr p N (W ntt ) = j^-tu / p N (w t (x)) dx 

n " W J[0,27r)" 

for all N G N. This, however, follows immediately from Lemma [2. 1[ 

The assertion about the differentiability follows from (fl71) . and a straightforward com- 
putation yields 



v w 2^ J \ 1 - rifxj'^tx) 1 -' y 



which is strictly positive unless ri(x) = r 2 (x) for almost every x. □ 

Recall the definitions of the Chernoff and the Hoeffding distances and their mean versions 
in Section 12. 1[ The above Lemma yields the following: 

3.7 Proposition. The mean Chernoff distance exists, and 

If supp pai,yi < supp p^mi ^ G N, then the mean Hoeffding distances exist, and 

—l r — 

H M {r\ p aim || Pa 2 m) = SU P — i : — ' r > °- 

o<t<i 1 — t 

Proof. The assertions follow immediately from the uniform convergence established in Lemma 
EJ3 □ 
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For operators < A, B < I on a finite-dimensional Hilbert space /C, define 

S 2 (A\\ B) := A (log A - log B) + (I — A) (log (I — A) — log (/ - B)) , 

which is a formal generalization of the relative entropy 5*2 (a 1 1 b) of the Bernoulli distributions 
(a, 1 — a) and (6, 1 — 6) that we get when dim/C = 1. 



(18) 
(19) 



3.8 Lemma. Assume that q\ and q 2 are strictly positive. Then, 

limi^-^l) = ^ Al + <?i(x)) S 2 (ri(x) || r 2 (x))dx = cT^(l) 

1 1 /" 27r 

lim — <9+</v(0) = -7T (l + ? 2 (x))S 2 (r 2 (x) || n(x))dx = 9+^(0). 

Proo/. Define B n (t) := f t (A^) + f x _ t (A 2 n) ) and /i n (t) := logc n (t) = ^(B^t^y, y). By 
a straightforward computation, 

h' n (t) = 2n(B n (ty 1 B' n (t)B n (ty 1 y,y) = AK{B n {t)~ l [X 1>n>t log R l>n - X 2 , n ,i_ t log # 2 , n ] ^(t)" 1 ^, 
where 

x M>t := (4 n) + - 4)7 [(4 n) + - (4 n) - 4)f , 

and one can easily see that 

d~h n (l) = \imti n (t) = -k (log R 2 , n y, y) ■ 

t/1 

A somewhat lengthy but otherwise again straightforward computation yields 
^ logdet ((/„ - W n , t ) = Tr [\ogR hn - log R 2 , n ] [l n - (I n - W n , t y l ] , 

which, in the limit t / 1, yields 

d~ log det (4 - W n , t ) | = Tr Q< n) [log R %n - log fl 1)Tl ] . 



Thus, by (1T4]). 
9>n(l) = W B (t) 

c/ 1 



- Tr log (/„ + gS n) ) + Tr log (l n + Q 2 n) ) - Tr [log R 2 , n - log i? 1>n ] 
-« (log R 2>n y, y) 



TrQ^ loggr-log^' 



(n) 



Tr (/ n + g[ n) ) [log (/ n + gS n) ) - log (4 + Q [ 2 n) ) 

-k {log R 2 , n y,y) (20) 
Tr (4 + Q i i ) )S 2 (R hn \\R 2>n ) - k (log R 2 , n y, y) ■ 
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Since the sequence log-R 2jn , n G N, is bounded, lim n h n {l) = 0, and hence, Lemma EOl 
applied to (I2P1) yields 

1 1 f 27T 

lim— d Vn(l) = tt - / <?i(x) [log gi (x) -logg 2 (x) 

- (1 + 9l (x)) [log(l + gi(x)) - log(l + g 2 (x))] rfx 
7r (l + ?i(x))5 2 (r 1 (x)||r 2 (x))dx 



= 0^(1), 

where the last identity follows by a straightforward computation from (1171) . This proves 
(ITS]) , and (1191) follows by a completely similar computation. □ 

3.9 Proposition. Assume that q\ and g 2 are strictly positive. Then, the mean rela- 
tive entropies SM(p ai , yi \\ Pa 2 ,y 2 ) an d SM(p a2 , y2 \ \ Pauyi) an d th e mean Hoeffding distance 
Pctim II Pa 2 ,y 2 ) exist, and 

1 f 2n 

SM(p ai , yi \\p a2 ,y 2 ) = —J (1 + 9i(x)) S 2 (ri(x) ||r 2 (x))dx = 9 ^(1) = H M (0\p auyi \\p a2 , y2 ), 

1 /" 27r 

5 W(P«a, w ||Pai, w ) = 7^ J (1 + g 2 (x)) S 2 (r 2 (x) 1 1 r^x)) dx = -<9 + ^>(0). 



Proof. The assertions follow immediately from Lemma l3T8l by the identities S(p$ >yi \ \ p^,% 
d~Ml) = H(0\ pg yi 1 1 pS, 2 ) and S(p^ 2 \ y2 1 1 pg yi ) = -<9+^„(0). □ 



,V2, 



3.2 Error exponents 

By the previous section, ip(t) = \im. n \ , ^ n {t) exists on [0,1]. Being the limit of convex 
functions, i/i is convex as well, and the uniformity of the convergence ensures that ip is also 
continuous on [0, 1]. Moreover, ip(t) < 0, t G [0, 1]. The polar function of ip is 

tp(a) := sup{ta - ij}(t) : t G [0, 1]}. 

For each a G IR and n G N, define the functions 

e n , a (T) := e-^ a a n (T)+/3 n (T) = e"^-Tr [e'^p^ - pgj T, T e T (B (F(H n ))) . 

Here we use the notation 

T{C) := {T G C : < T < 1} 

to denote the set of tests in a C*-algebra C. Let S^a := {e~ n " a pai,y — Pa 2 ,y > 0} be 
the spectral projection corresponding to the positive part of the spectrum of the self-adjoint 
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operator e~ n " a pa 1 \ y — pa 2 ,y- S n ,a is usually referred to as the Neyman-Pearson test or Holevo- 
Helstrom test, and is easily seen to be a minimizer of e n _ a on T (B (jF(7Y n ))). 

Theorem 1 in [1] tells that for positive semidefinite operators A and B on some Hilbert 
space /C, 

- Tr(A + B) - - Tr \A - B\ < Tr A^ 1 " 1 , t G [0, 1]. 

Substituting A := e~ nVa p^ y and S := /3^ y , we get 

e n , a (Sn, a ) < e~ tn " a Tr (p^J (pSJ 1 "*, * G [0,1], 
and hence, 

limsup — loge^S^a) < inf {-to + ^(t)} = -cp(a). (21) 

n U u 0<i<l 

Note that SVi.a is not necessarily in the observable algebra CCR(7Y n ) of the local system 
on the cube C n . However, we have the following: 

3.10 Lemma. For all a G R, there exists a sequence of tests S nja G T(CCR(7i n )), n G N 
such that 

limsup \e nA {S nA ) < -f{a), (22) 

n— >oo 71 

limsup — log an(^n, a ) < ~{v( a ) ~ «}> (23) 

n— >oo 71 

limsup^ log p n (S nta ) < -<p(a). (24) 

n— »oo 71 

Proof. We use that the von Neumann algebra generated by CCR(7Y n ) is equal to B (J r (TC n )) 
[5], Proposition 5.2.4], and hence it contains S n>a for all a G R. Then, by Kaplansky's density 
theorem, there exist tests S nsa G T(CCR(7Y n )) such that 

Hence, (1221) follows by ( 12~T|) . and the rest is immediate from (5{S n>a ) < e nja {S n>a ) and 
e~ n a a n (S n 



For the rest, we will rely on the analysis in the paper [13] . Note that our setting here 
is somewhat different from that of [13], as the local algebras are infinite-dimensional and 
the scaling in the asymptotics is \ instead of -. However, most of the analysis in [T5] 
carries through whenever the existence and differentiability of x/i on [0, 1] can be established. 
Below we show how the results of [13] can be adapted to the present setting. Note that 
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there are different sign conventions in the literature in defining the error exponents. Here 
we chose a sign convention opposite to the one in [T3] and used a different notation in order 
to emphasize the equality of the error exponents and the corresponding relative entropy-like 
quantities. The correspondence between the notations of this paper and those of [13] is 

h{r\p ai , yi \\pa 2 , y2 ) = -B_(r\ p\\ a) , s {p ai , yi \ \ Pa 2 ,y 2 ) = -B(p\\a) , 
h{r\ p auyi \\pa 2 , y2 ) = -B(r\p\\a), s(p aim \\ p a2tV2 ) = -B(p\\a), 
h ( r \ P ai , yi \\ Pa 2 , y2 ) = -B(r\p\\a) , s(p auVl \\p a2jV2 ) = -B(p\\a) , 

with v = 1 omitted from the notations. 

3.11 Lemma. Assume that qi and q 2 are strictly positive. For any sequence of tests T n e 
T(CCR(H n )), n E N, and d + ip(0) < a < d~ip(l), 

liminf —e na (T n ) > -ip(a). 

rwoo fl 

Proof. As S n>a minimizes e n ^ a over T (£> (JF(7Y„))), it is enough to show that 
liminf —e na (S n , a ) > -<p(a). 

n— >oo 72 

Let 

Ai |B (m, m ) := — log , X 2 , n (m, m):=— log - 

n p ljn (m, m'J n u p 2 ,n [ZH-, HE ) 

be random variables on J n . As it was shown in [19], 

2e n , Q (S n , a ) > e- n "> lin ({X 1>n > -a}) + p 2 , n {{X %n > a}) , (25) 

(see also the proof of [T3l Theorem 3.1]). The logarithmic moment generating function of 
X n>k with respect to pk, n is 

$ M (t) := log f e tx d{p Kn o A"-*) = log ^ e ^(s,m') pk n ^ ? 

R (m,m')eJ„ 

and one can immediately see that 

$i, n (^) = log p ljn (m^m ) 1 ^P2,n (m,2l')* = ^( 1 _ *)> 

(m,m')EJ„ 

$2,n(^0 = log ^ pi^^mf P2,n(Ek,m ) 1 ~ t = ^Pn(t)- 

(m,rn')£j n 
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By Lemma IBToT the asymptotic logarithmic moment generating functions $&(£) := lim n \<&k,n 
exist on [0, 1] and are differentiable in (0, 1), with & 1 (t) = ip{l - t), $ 2 (*) = fp(t), t G [0, 1]. 
The Gartner- Ellis theorem (see, e.g. [8]) then yields 

liminf — pi, n ({Xi,n > ~ a }) — ~ SU P { — at — <3?i (t)} = — sup {— at — — t)} = a — <p(a) 
n n u o<t<i o<t<i 

and 

liminf — p 2 ,n ({AT 2 ,n > a}) > - sup {at - $2 CO} = - sup {at - ip(t)} = —<p(a) 
n n u o<t<i o<t<i 

for all d + ip{Q) < a < d~ip{l). Thus, by (1251) . the assertion follows. □ 

3.12 Corollary. Assume that q± and q 2 are strictly positive. For any d + ip(0) < a < d~ip(l), 
there exists a sequence of tests S Uja G T(CCR(TC n )), n G N, such that 

lim — e na (S na ) = -<p(a), (26) 
n n u 

lim — loga n (£ re)0 ) = -{<p(a) - a}, (27) 

n n u 

lim — log/3 n (S n)0 ) = -cp(a). (28) 
n n v 



Proof. The first assertion follows immediately from Lemmas l3.10l and l3.lll Lemma 4.4 in [13] 
tells that lim inf n ^7 log (5 n (T n ) > — <p(a) for any sequence of tests for which lim sup n ^7 log a n {T n 
< —{if (a) — a} holds, and hence (I28I) follows from (123]) . Similarly, (1271) follows from (1241) by 
Remark 4.6 in [13]. □ 



Now we are in a position to prove our main result: 
3.13 Theorem. Assume that q\ and q 2 are strictly positive and q\ ^ q 2 . Then, 

c (.Pai,yi II Pai,y-z) c(Pqi,j/i || Pa 2 ,y 2 ) c (Pai,j/i 1 1 Pa 2 ,y 2 ) ^M\Pa\,yi 1 1 Pa 2 , y2 )i (29) 

M r |Pai,J/i ||Pq!2,2/2) = — ( r I Pai,yi \ \ Pai2,y2) = h(r\ p a \,y\ \\Pct 2 , y2 ) = -^Af ( r I Paa ,J/i II Pa.2,Vz)i (30) 
s (Pai,j/i || Pa 2 ,y 2 ) = §.(Pa 1 ,y 1 \ \ Pa 2 ,y 2 ) = s (Pa 1 ,y 1 \ \ Pa 2 ,y 2 ) = (Pai || Pa 2 ,y 2 )i (31) 

where ( 130]) holds for all < r < SM(p a2 , y2 II P<*i,3/i)- 

Proof. The assumptions yield that ^(0) = "0(1) = 1 an d <9 +, ?/>(0) < < 9 _, 0(1). Hence, by 
choosing a = in Lemma 13. 11] we get 

C(p auyi \\pa 2 ,y 2 ) < <f(0) = -mm#) = C M (p aim \\p a2 ,y 2 ). 
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On the other hand, Corollary 13.121 ensures the existence of a sequence of tests 5„ )0 , n G N 
such that lim n ^ 7 e rii o(5'n,o) — — <p(fy, an d hence, c(p ai>yi \ \ p a2 ,y 2 ) — fifyj from which f[2"9"j) 
follows. 

To prove (1301) . define </5>(a) := y?(a) — a, a G R. By Lemma 4.1 (see also Figure 2) in [T3j . 
tp is strictly monotonically increasing on [d + ip(0),d~ip(l)}, with range [0, d~ip(l)], while </? 
is strictly monotonically decreasing on the same interval, with range [0, —d + tp(0)]. Hence, 
for any < r < — d + tp(0) = Sjw(/? a2)Wa II Pai,yi)> one can find a unique a r G (9 +, 0(O), 
such that <p(a r ) = r. By Corollary 13.121 we have for any d + ip(0) < a < a r , 

lim — logo: n (£ re)0 ) = -0(a) < -<p(a r ) = -r and lim — log/3 n (£ n;(1 ) = -<p{a). 
n n n n 

Hence, 

K T \ Paim II Pa 2 ,y 2 ) > SUp (f(a) = <f(a r ), 
d+i}(0)<a<a r 

and the rest of the proof goes exactly the same way as in the proof of Theorem 4.8 in [T3] . 

The last assertion follows immediately from Propositions 5.1 and 5.2 in [13], by noting 
that in our setting, ip(l) = and d~^{l) = Sm^oh^ II Pai,yi)- ^ 



4 Conclusion 

We considered the hypothesis testing problem of discriminating two Gaussian states of an 
infinite bosonic lattice, and gave complete solutions to the problems of the Chernoff bound, 
the Hoeffding bound and Stein's lemma under the assumptions that both states are gauge- 
invariant with translation-invariant quasi-free parts. 

Note that the natural structure underlying the theory of the CCR algebra and Gaussian 
states is a real vector space H, equipped with a symplectic form. On the other hand, if H 
is finite dimensional then a Gaussian state always defines a canonical complexification of H 
in which the state becomes gauge- invariant. Our assumption that both states are gauge- 
invariant can heuristically be understood as requiring that the two states yield the same 
complexification, which is clearly the strongest technical limitation of our approach. It is an 
open question how to extend our results to the state discrimination problem for non-gauge 
invariant Gaussian states. 

To be able to treat the infinitely extended lattice, we have chosen a C*-algebraic descrip- 
tion of the system. In particular, we defined states of the system as linear functionals on the 
observable algebra, as a density operator may not exist in this case. In all computations, 
however, we used a concrete representation of the CCR algebra, the Fock representation. 
Lemma I'd. 101 shows that, as far as our asymptotic state discrimination problem is concerned, 
it does not matter what representation we use, as the asymptotically optimal performance 
can be reached by measurement operators from the CCR algebra. 
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Appendix A 

If A G B(H.) then A® m leaves V m 7i invariant, and we denote its restriction to \J m 7i by V m A 
The Fock operator A F , corresponding to A, is 



oo 



A F : = V m A with V(A F ) := I ®™ =0 x m G F{H) : £ ||(V"M)x m f < oo 

m=0 I. m=0 

Note that the Fock operators are closed, and 

T S {H) := {®% =0 x m : x m G y m H < m < M, M G N} 

is a common core for all Fock operators, on which ApBp = (AB)p holds. If A > then we 
also have (Ap) 1 = (^4*)^ on jFy(7i) for any t G R, with the convention 0* := 0, t G R. Fock 
operators are also characterized by the property A F xp = (Ax) F , x £TC. 

If A > is a finite-rank operator and A = ^I-=i ^fcl e fc)(efc| is an eigen-decomposition of 
A, then 

V m A= ^ A„|e„)(e M | (32) 

m\+ . . .+m r =m 

is an eigen-decomposition of \/ m A, where 

X m := XT 1 ■ . . . ■ Xr, := - 1 V f/( m ) e f mi ® . . . ® e f*-, (33) 

\/m\\ . . . m r \m\ *—r 

and Ua , cr G S^, denotes the standard unitary representation of the symmetric group S m 
on H® 171 . As a consequence, 

oo r / oo \ 

£TrV m A = n PTA™ , (34) 

m=0 fc=l \m=0 / 

which is finite if and only if A < I, in which case Ap is trace-class with 

TrA F = det (I -A)' 1 . (35) 

If T(B) := ©^ =o r m (_B) is the second-quantized version of a finite-rank operator B G B{7i), 
where Y m {B) is the restriction of Y!k=i I® {k ~ l) ®B® I m ~ k onto V m ft, then 

Tr ™ = dS(7r7) Tr 7^4 B - (36 » 
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